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G\ ! ABSTRACT 

m . 

Context. The star-to-star scatter in lithium abundances observed among otherwise similar stars in the solar-age open cluster M 67 is one of the 
most puzzling results in the context of the so called "lithium problem". Among other explanations, the hypothesis has been proposed that the 
dispersion in Li is due to star-to-star differences in Fe or other element abundances which are predicted to affect Li depletion. 
Aims. The primary goal of this study is the determination of the metallicity ([Fe/H]), a- and Fe-peak abundances in a sample of Li-poor and 
Li-rich stars belonging to M 67, in order to test this hypothesis. By comparison with previous studies, the present investigation also allows 
J—!, us to check for intrinsic differences in the abundances of evolved and unevolved cluster stars and to draw more secure conclusions on the 
Or abundance pattern of this cluster. 

^ . Methods. We have carried out an analysis of high resolution UVES/VLT spectra of eight unevolved and two slightly evolved cluster members 
'_h ' using MOOG and measured equivalent widths. For all the stars we have determined [Fe/H] and element abundances for O, Na, Mg, Al, Si, Ca, 
^ \ Ti, Cr and Ni. 

■ Results. We find an average metallicity [Fe/H] = 0.03 ± 0.01, in very good agreement with previous determinations. All the [X/Fe] abundance 
ratios are very close to solar. The star-to-star scatter in [Fe/H] and [X/Fe] ratios for all elements, including oxygen, is lower than 0.05 dex, 

• i— l , implying that the large dispersion in lithium among cluster stars is not due to differences in these element abundances. We also find that, when 

■ using a homogeneous scale, the abundance pattern of unevolved stars in our sample is very similar to that of evolved stars, suggesting that, at 
?H ' least in this cluster, RGB and clump stars have not undergone any chemical processing. Finally, our results show that M 67 has a chemical 

. . . i composition that is representative of the solar neighborhood. 



Key words. Stars: abundances - Stars: Evolution - Open Clusters and Associations: Individual: M 67 



1. Introduction 

"Classical" or "standard" models of stellar evolution include 
convection only as internal mixing process and neglect more 
complex physical processes like diffusion, magnetic fields, ro- 
tation; these models predict that solar-type stars should not de- 
plete lithium (Li) while on the main sequence (MS), since the 
base of the convective zone does not reach the layers in the 
stellar interior where the temperature is high enough for Li re- 
actions. Furthermore, according to standard models, stars with 
the same age, mass, and chemical composition should undergo 
the same amount of Li depletion. In sharp contrast with these 
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predictions, not only is now well established on observational 
grounds that solar-type stars, including the Sun, do deplete a 
significant amount of Li during the MS phases, but Li deple- 
tion can be different for similar stars (e.g., Randich 2006 and 
references therein). The factor of about 10 star-to-star scatter 
in Li abundances observed among otherwise identical F- and 
G-type members of the 4.5 Gyr old cluster M 67 indeed rep- 
resents one of the most puzzling results in the context of the 
so-called "Lithium problem" (Spite et al. 1987 Garcia Lopez 
et al. 119551 Pasquini et al. H9971 Jones et al. [T999l . 

A large dispersion in Li abundances is present among old 
field stars (e.g., Pallavicini et al. 119871 Pasquini et al. 1994 1. On 
the other hand, whereas no dispersion is seen in the 600 Myr 
old Hyades, in the three 2 Gyr old clusters IC 4651, NGC 3680 
and NGC 752 (Randich et al. l2TjQTjl Sestito et al.|T004), nor 
in 6 Gyr old NGC 188 (Randich et al l2003l . preliminary re- 
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suits of Li measurements among large samples of stars in the 
2 Gyr, metal rich NGC 6253 and in the very old Collinder 261 
suggest that they could be characterized by some amount of 
scatter (Pallavicini et al. 120061 Randich |2006> . In other words, 
the appearance of a dispersion seems to depend more on the 
characteristics of the cluster than on age. 

It has been shown that non membership and/or binarity are 
not the reasons for the scatter in M 67 (Pasquini et al. [1997 1; 
also, the effects of chromospheric activity on the formation of 
the Li i line, which are proposed as a possible explanation for 
the dispersion in young clusters (Jeffries 2006 and references 
therein), are unlikely to play a role, given that at the old age of 
M 67 the level of chromospheric activity should be low enough 
not to affect the Li i line. Hence, the scatter is most likely intrin- 
sic and, under the very reasonable assumption that cluster stars 
were all born with the same Li content, it must reflect different 
amounts of Li depletion. 

Different possibilities were proposed to explain the exis- 
tence of the spread (Randich 2006 1 ; among them, it has been 
suggested that M 67 members do not have an identical chemi- 
cal composition and that the scatter in Li may reflect a scatter in 
iron content or, more in general, in heavy element composition 
(Garcia Lopez et al. fT988l Piau et al. l2"00"3l . 

As well known, chemical composition affects stellar opaci- 
ties and thus internal structure, mixing processes (both standard 
and non standard ones) and, in principle, Li depletion. The ef- 
fect of variations of the chemical composition on pre-main se- 
quence Li depletion has been theoretically investigated in dif- 
ferent studies (e.g, Swenson et al. 1 19921 Piau & Turck-Chieze 
2002, Sestito et al. 2006 1 and all of them agree in that even 
relatively small changes of the mass fraction of elements crit- 
ical for the opacity can have significant effects on the amount 
of Li depletion. Similarly, Piau et al. (2003 1 have shown that 
variations in CNO abundances can change the amount of Li 
depletion during the MS phases and suggested that star-to-star 
differences in these element abundances or, more in general in 
a and Fe-peak abundances, could explain the observed scatter 
in M 67. More quantitatively, they found that a difference of 
0.05 dex in [CNO/Fe] would result in a difference in log n(Li) 
of ~ 0.5 dex (i.e, smaller than the observed spread), implying 
that larger differences in CNO abundances are needed to ex- 
plain the whole spread. 

Garcia Lopez et al. ( 1988 1 did not find any difference in the 
overall metallicity of Li-poor and Li-rich stars and ruled out 
that different [Fe/H] values within M 67 could be the reason 
for the scatter in Li. Here we extend their work to elements 
other than iron and to a larger sample of stars, to investigate 
a) whether a large (> 0.05 dex) star-to-star scatter in heavy 
elements is present among unevolved cluster members; and, in 
case such a spread is detected, b) whether it is related to the 
dispersion in Li. 

In a more general context, the determination of the metal- 
licity and chemical composition of open clusters covering a 
large interval of ages, metallicities, and Galactocentric dis- 
tances is a critical tool to investigate the formation and evolu- 
tion of the Galactic disk (Friel 2006 and references therein). 
M 67 is one of the closest and best studied old open clus- 
ters; nevertheless, relatively few studies focused on the deter- 



mination of its chemical composition (Tautvaisiene et al. 2000 
and references therein), and most of them are based on few 
and/or evolved stars that may have undergone chemical pro- 
cessing. The present work represents the first abundance study 
of M 67 based on a rather large sample of unevolved (or slightly 
evolved) cluster members, allowing us to check whether in- 
trinsic differences exist between abundances of unevolved and 
evolved stars and to draw more secure conclusions on the abun- 
dance pattern of this cluster. The paper is structured as follows: 
in Sect. 2 the sample and the observations are described, while 
the abundance analysis is presented in Sect. 3. The results, to- 
gether with a discussion of internal and systematic errors, and 
a comparison with findings from previous studies are presented 
in Sect. 4. Finally, a discussion of the results and conclusions 
are given in Sect. 5. 



2. Sample and observations 

Our sample includes 10 single M 67 members that are listed 
in Tabled together with V, B - V, and (B - V)o colors. The 
color-magnitude diagram of the cluster with our sample stars 
evidenced in black is shown in Fig. the diagram indicates 
that seven of them are still on the MS, one is close to the the 
turn-off, while the remaining two are already evolved to the 
subgiant branch. 

The observations were obtained using the UVES spectro- 
graph (Dekker et al. 1200 2 1 on VLT UT2/Kueyen during three 
different observing runs; the first one was carried out in Visitor 
mode in April 2000, while the other two were performed 
in Service mode in Fall 2001 and Spring 2002, respectively. 
Details on the observations, whose primary goal was the mea- 
surement of beryllium abundances in the sample stars, were 
already given in Randich et al. (2002 1 and we summarize the 
main points in the following. UVES was operated in Dichroic 
Mode using Cross Dispersers #1 and #3 in the Blue and Red 
arms, respectively. The Red arm, which is of interest here, is 
equipped with a mosaic of two CCDs composed by a EEV 
2048x4102 CCD and a MIT-LL 2048 x 4102 CCD; the spec- 
tral coverage ranges from approximately 4780 to 6810 A. The 
15/j.m pixels together with a 1 arcsec wide slit (projecting onto 
4 pixels) and CCD binning lxl, yielded a resolving power 
R~ 45, 000. Exposure times were set based on the requirement 
of a good S/N ratio in the near-UV spectral regions where the 
Be n lines are located and range between 1.3 and 3 hrs per star. 
Data reduction was carried out using the UVES pipeline and 
following the usual steps. Typical S /N ratios per resolution el- 
ement measured on the extracted 1-D spectra range between 90 
and 180. 



3. Abundance analysis 

Abundance analysis was carried out by means of measured 
equivalent widths (EWs) and using Version 2000 of MOOG 
(Sneden 1973 1 with a grid of 1-D model atmospheres from 
Kurucz A1993I . We recall that MOOG performs a standard-LTE 
analysis. 
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Table 1. Sample stars, photometry, stellar parameters and derived metallicities together with rms scatter. Numbering is from 
Sanders ( 1977 ». V magnitudes and B - V colors were taken from Montgomery et al. J1993I . (B - V) colors were instead 
retrieved from Jones et al. ( 1999 1 who had derived them based on original B - V and V - I colors of Montgomery et al. and 
assuming a reddening E(B - V) — 0.05. For stars S 1034 and S 1239, not included in the sample of Jones et al. ( 1999 1, we obtained 
(B - V)o colors in the same way. One Hyades member is also included (see text). 



S V B-V (B-V)o T eff £, log g [Fe/H]i rms 











(K) 


(km/s) 




969 


14.18 


0.665 


0.622 


5800 


1.10 


4.4 


988 


13.18 


0.570 


0.534 


6153 


1.45 


4.1 


994 


13.18 


0.581 


0.535 


6151 


1.45 


4.1 


995 


12.76 


0.559 


0.521 


6210 


1.50 


3.9 


998 


13.06 


0.567 


0.518 


6223 


1.50 


4.0 


1034 


12.65 


0.608 


0.567 


6019 


1.50 


4.0 


1239 


12.75 


0.758 


0.692 


5541 


1.25 


3.8 


1252 


14.07 


0.643 


0.588 


5938 


1.15 


4.4 


1256 


13.67 


0.660 


0.595 


5907 


1.10 


4.5 


2205 


13.14 


0.566 


0.534 


6156 


1.45 


4.1 


vB 187 


9.01 


0.76 


0.75 


5339 


0.9 


4.5 



3.1. Line list and equivalent widths 

Spectral lines to be used for the analysis of Fe i, Fe n and, other 
elements (Na, Mg, Al, Si, Ca, Ti, Cr, Ni) were selected from 
different sources in the literature and subsequently checked for 
suitability (in particular for blends) on the solar spectrum ob- 
tained with UVES at the same resolution as our sample stars. 
We finally retained in the list 55 Fe i and 1 1 Fe n lines, and a 
total of approximately 70 lines for all the other elements. 

The task eq in SPECTRE was used to interactively measure 
the EWs of the spectral lines by gaussian fitting. Although the 
spectra had been previously normalized, local continuum was 
inspected and, if needed, adjusted at each EW measurement. 

3.2. Atomic parameters and solar analysis 

The majority of log gf values for both Fe and other element 
lines were taken from the Vienna Atomic Line Data-base 
(VALD: Piskunov et al. [T995 :; Kupka et al 1999 Ryabchikova 
et al. 119991 — http://www.astro.uu.se/htbin/valdi. For a few 
lines for which log gf from VALD were not available or gave 
very discrepant abundances for the Sun, log gf were either 
taken from other sources in the literature or, in a few cases, ad- 
justed through an inverse solar analysis. Note that log gf values 
for all Fe i lines were taken from VALD. In Table 2 we show all 
the lines included in our list together with their log gf values, 
the source for these values, and the EWs measured in the solar 
spectrum. Radiative and Stark broadening are treated in a stan- 
dard way in MOOG; as for collisional broadening, we used the 
Unsold approximation ( 1955 1 for all the lines. As discussed by 
Paulson et al. (2003) this choice should not greatly affect the 
differential analysis with respect to the Sun. We also mention 
that very strong lines that are most affected by the treatment of 
damping have been excluded from our analysis. 

The analysis of the solar spectrum was performed assuming 
the following solar parameters: T e ff = 5770 K, logg Q = 4.44 



0.01 


0.04 


0.03 


0.04 


0.0 


0.04 


0.05 


0.03 


0.07 


0.05 


0.01 


0.05 


0.02 


0.03 


0.05 


0.03 


0.06 


0.03 


0.00 


0.05 


0.13 


0.05 



and £ = 1.1 km/s. We mention that log n(Fe i) vs. EW did 
not show any trend, implying that the assumed microturbulence 
is correct. On the contrary we found a small, but statistically 
significant (slope equal to 0.009, correlation coefficient 0.326) 
trend of Fe i abundances vs. excitation potential (EP). Such 
a trend would disappear by assuming a 70 K higher solar ef- 
fective temperature (T e ffo = 5840 K), which would result in 
larger and smaller abundances for Fe i and Fe n, respectively 
(log n(Fe i)=7.58 and log n(Fe n)=7.49). 

Output solar abundances are listed in Table[3]together with 
those from Anders & Grevesse ( 1989) that are used as input in 
MOOG. Note that, although solar abundances from Anders & 
Grevesse ( 1989 ) have been superseded by the study of Grevesse 
& Sauval ( 1998), the two studies yield the same abundances or 
very little differences for the elements analyzed in this study. 
Table[3]shows a good agreement between the abundances deter- 
mined by us and the values of Anders & Grevesse ( 1989), the 
most discrepant element being Na with Alog(Na)=0.05 dex. 
We also find a very good agreement between log n(Fe) from 
Fe i and Fe n. 



3.3. Oxygen 

We do not discuss here all the intricacies related to oxygen 
abundance determinations in stars, and we refer to Bensby et al. 
(2004) and Schuler et al. (2005 1 for recent discussions. The for- 
bidden lines at 6300.30 A and 6363.78 A are not significantly 
sensitive to NLTE effects nor to stellar effective temperature; 
therefore, as widely acknowledged, these lines are the most re- 
liable ones for deriving O abundances. Our O determinations 
are based on the [O i] forbidden line at 6300.3 A. Although 
this line is located in a spectral region that is severely affected 
by the presence of telluric lines, the radial velocity of M 67 is 
such that they do not affect the [O i] line in our spectra. 
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Table 3. Solar abundances from Anders & Grevesse ( 1989 1 and 
those determined through our analysis. Note that logg/ values 
for Al were derived by an inverse solar analysis. 



Element 


log n(X) A G89 


log n(X) our 


Oi 


8.87 


8.66±0.04 


Nai 


6.33 


6.28±0.03 


Mgi 


7.58 


7.57±0.01 


All 


6.47 


6.47±0.03 


Si i 


7.55 


7.56±0.03 


Cai 


6.36 


6.35±0.02 


Til 


4.99 


4.97 ±0.02 


Cri 


5.67 


5.65±0.02 


Fei 


7.52 


7.52±0.03 


Fe ii 


7.52 


7.51 ±0.04 


Nil 


6.25 


6.25 ±0.03 



As first pointed out by Lambert dl978> . the 6300.3 A line 
is blended with a Ni i feature whose contribution must be taken 
properly into account for an accurate determination of O abun- 
dance. Hence, we determined O both for the Sun and our sam- 
ple stars using the driver blend in MOOG, that allows account- 
ing for the blending Ni i 6300.34 feature when force fitting the 
abundance of the O to the measured EWs of the 6300.3 A fea- 
ture. As done for the other lines, EWs of the 6300.3 A feature 
were measured by gaussian fitting using SPECTRE. The solar 
[O i] EW measured in the UVES spectrum is 5.5 +0.3 mA, 
in good agreement with previous estimates (see Schuler et al 
2005 1; EWs for the sample stars are listed in the second col- 
umn of Table |5] 

For our analysis we employed the [O i] g/-value deter- 
mined by Allende Prieto et al. (1200 1> . log gf = -9.717, while 
for the Ni 1 6300.34 A blend, following Johansson et al. ( 1217031 . 
we used log gf — -2.1 1. Note that, as discussed by Schuler et 
al., we would have obtained virtually the same O abundances 
by considering the two separate components of the Ni i blend. 
Assuming a solar Ni abundance log n(Ni)=6.25 (see Tabled, 
we obtained for the Sun log n(0)=8.66 +0.04, much below the 
classical value of Anders & Grevesse ( 1989 1, but in good agree- 
ment with recent determinations, including those obtained us- 
ing 3-D analysis (see again the discussion in Schuler et al.). 

3.4. Stellar parameters 

Initial stellar parameters for the sample stars were estimated 
from photometry (see Table [Q. Effective temperatures were 
retrieved from Jones et al. (I1999> . who, in turn, had derived 
them using the calibration of Soderblom et al. (1993 1, i.e, 
T eff = 1 808 x (B - V)2 - 6 1 03 x (B - V) + 8899. Microturbulence 
velocities were estimated as: £ = 3.2 x 10" 4 (T eff - 6390) - 
1.3(logg - 4.16) + 1.7 (Nissen 1981), while surface grav- 
ities were derived using the relationship: logg = 4.44 + 
logM/M + 41ogT eff /T eff + 0.4(M bol - M bol 0). Absolute 
bolometric magnitudes were determined from V magnitudes, 
adopting a distance modulus (m-M)o=9.6 (Pace et al. [2004 
Sandquist I2004> and assuming bolometric corrections from 



Johnson ( 1966 1. Finally, stellar masses were derived following 
Balachandran dl995> . 

In order to determine spectroscopic temperatures one 
would need to change the initial T e ff values until no log n(Fe)i 
vs. EP trend is seen. However, we decided not to change ini- 
tial T e ft derived from photometry, since a trend is present for 
the solar spectrum and for all the sample stars we found simi- 
lar slopes (in magnitude and sign). To further test the correct- 
ness of our approach, for each star we plotted Alog n(Fe)= 
log n(Fe) -log n(Fe) stal as a function of EP and checked that no 
significant residual trends were present. Spectroscopic micro- 
turbulence values for each star were instead obtained by chang- 
ing the input microturbulence until no log n(Fe) vs. EW trends 
were seen. Finally, log g values were varied only when nec- 
essary, in order to have a difference between iron abundances 
derived from Fe i and Fe n lines not greater than 0.05 dex, 
i.e., consistent, within errors, with the ionization equilibrium. It 
turned out that we needed to change the initial values of gravity 
for two stars only, S1239 and S1256. Assumed parameters for 
the sample stars are listed in Cols. 5-7 in Tabled 

4. Results 

4.1. Abundances 

Log n(X) values for each line were determined based on mea- 
sured EWs and stellar parameters listed in Tabled Final abun- 
dances for each star and each element were determined as the 
mean abundance from the different lines, lcr clipping was per- 
formed for iron, but not for the other elements for which a 
smaller number of lines was available. As to oxygen, [O/H] 
and [O/Fe] values for the sample stars were obtained from the 
measured EW of the 6300.3 A blend (see Table |5} and assum- 
ing Ni abundances determined through our analysis. As already 
mentioned, [Fe/H] and [X/Fe] ratios for each star were deter- 
mined differentially with respect to solar abundances listed in 
Tabled 

[Fe/H]i abundances and rms values for the sample stars are 
listed in Cols. 8 and 9 of Table ^ while in Fig. [2] we show 
[Fe/H] as a function of effective temperature. The figure clearly 
shows that no trend of [Fe/H] as a function of T e g is present. 
Using all the stars in our sample, we obtain a mean value 
[Fe/H]= 0.03 + 0.01 (rms= 0.03, 10 stars), very close to so- 
lar. This mean value is listed in the first column of Table0 

The [X/Fe] ratios for Na, Mg, Al, Si, Ca, Ti, Cr and Ni 
are listed in Table0] Errors on [X/Fe] values in Table 0] corre- 
spond to the quadratic sum of rms of [Fe/H] and [X/H] values. 
Results for O are listed separately in Table|5] where both [O/H] 
and [O/Fe] values are given; in this case errors in [O/H] corre- 
spond to uncertainties in the measured EWs of the forbidden 
line, while errors in [O/Fe] are the quadratic sum of errors in 
[O/H] and rms of [Fe/H]. Mean abundance ratios for M 67 to- 
gether with 1 cr standard deviation are listed in the second col- 
umn of Table [7] while [X/Fe] ratios for all elements are plot- 
ted in Fig. [3] as a function of effective temperature. The figure 
shows that, as in the case of [Fe/H], no evident trends of [X/Fe] 
ratios with T e ff are present. A small amount of star-to-star scat- 
ter might be present, both in [Fe/H] and [X/Fe] ratios, but for 
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Table 4. Derived abundance ratios for M 67 stars and one Hyades member. For each star, quoted errors on abundance ratios are 
the quadratic sum of rms values of [Fe/H] and [X/H]. Average abundance ratios for the Hyades, taken from Friel (2006|l, are 
given in the last line. 



s 


[Na/Fe] 


[Mg/Fe] 


[Al/Fe] 


[Si/Fe] 


[Ca/Fe] 


[Ti/Fe] 


[Cr/Fe] 


[Ni/Fe] 


969 


0.00 ± 0.06 


-0.01 ± 0.05 


-0.01 ± 0.07 


-0.01 ±0.06 


0.09 ± 0.06 


0.00 ± 0.07 


0.00 ±0.06 


-0.04 ± 0.07 


988 


0.12 ±0.04 


0.03 ±0.13 


-0.04 ± 0.06 


0.10 ±0.05 


0.12 ±0.06 


0.00 ± 0.06 


0.09 ± 0.04 


0.03 ± 0.06 


994 


0.12 ±0.05 


0.02 ± 0.09 


-0.04 ± 0.04 


0.02 ± 0.07 


0.07 ± 0.06 


0.04 ± 0.06 


0.03 ± 0.05 


0.04 ± 0.06 


995 


0.04 ± 0.04 


0.00 ±0.1 3 


-0.06 ± 0.08 


0.03 ± 0.06 


0.05 ± 0.05 


-0.06 ± 0.08 


-0.04 ± 0.04 


-0.07 ± 0.04 


998 


0.08 ± 0.06 


0.00 ± 0.09 


-0.08 ± 0.07 


-0.02 ± 0.07 


0.01 ± 0.07 


0.00 ± 0.08 


-0.05 ± 0.05 


-0.03 ± 0.06 


1034 


0.07 ± 0.06 


-0.02±0.05 


-0.08 ± 0.04 


0.03 ± 0.06 


0.03 ± 0.06 


0.00 ± 0.07 


0.00 ±0.05 


-0.03 ± 0.07 


1239 


0.03 ± 0.09 


-0.01±0.09 


0.02 ± 0.06 


0.05 ± 0.07 


0.00 ± 0.07 


-0.05 ± 0.05 


-0.01 ±0.05 


-0.03 ± 0.05 


1252 


-0.04 ± 0.06 


-0.04 ± 0.03 


-0.09 ± 0.06 


0.00 ± 0.05 


0.02 ± 0.04 


-0.08 ± 0.06 


-0.04 ±0.04 


-0.05 ± 0.05 


1256 


-0.06 ± 0.04 


-0.04 ± 0.04 


-0.11 ±0.05 


-0.03 ± 0.05 


0.04 ± 0.06 


-0.04 ± 0.05 


-0.02 ± 0.07 


-0.04 ± 0.05 


2205 


0.13 ±0.06 


0.03 ± 0.07 


-0.06 ± 0.08 


0.04 ± 0.08 


+0.03 ± 0.06 


-0.06 ± 0.09 


-0.01 ±0.05 


0.03 ± 0.06 


vB187 


+0.00 ± 0.05 


-0.06 ± 0.02 


-0.07 ± 0.07 


0.07 ±0.07 


0.04 ± 0.06 


-0.09 ± 0.04 


-0.05 ± 0.02 


0.03 ± 0.06 


Hyades li( 


0.01 


-0.06 


-0.05 


0.04 


+0.06 


-0.06 







Table 5. [O i] 6300 A line EWs, [O/H] and [O/Fe] values for the sample stars and the Hyades member. Quoted errors in [O/H] 
are due to uncertainties in EWs, while errors in [O/Fe] are the quadratic sums of errors in [O/H] and rms values of [Fe/H] (see 
Table^. In the last line we provide average [O/H] and [O/Fe] ratios for the Hyades from Schuler et al. (2005 1. 



s 


EW (6300.3 A) 




[O/H] 




[O/Fe] 




(mA) 










969 


5.6 J 


:0.5 


0.02 


±0.06 


0.01 


±0.07 


988 


5.3 + 


:0.9 


0.00 


+0.10 


-0.03 


±0.11 


994 


5.0 ± 


:0.9 


-0.03 


+ 0.10 


-0.03 


±0.11 


995 


6.6 ± 


:0.9 


+0.08 


±0.08 


0.03 


±0.08 


998 


5.7 4 


:0.5 


+0.04 


±0.04 


-0.03 


±0.06 


1034 


6.5 4 


:0.5 


+0.03 


±0.04 


+0.02 


±0.06 


1239 


10.2 


±1 


+0.03 


±0.06 


+0.01 


±0.07 


1252 


5.1 4 


:0.4 


+0.00 


±0.05 


-0.05 


±0.06 


1256 


5.0 4 


:0.5 


+0.03 


±0.05 


-0.03 


±0.07 


2205 


5.4 


±1 


+0.01 


±0.11 


0.01 


±0.12 


vB187 


5.6 4 


:0.4 


-0.19 


±0.08 


-0.32 


±0.09 


Hyadesschuic, 








0.19 




0.10 



all the elements the scatter is well within measurement uncer- 
tainties. A more detailed discussion on the possible presence of 
a scatter will be provided in Sect. 15.11 



4.2. Errors 

Sources of internal errors include uncertainties in atomic and 
stellar parameters, as well as errors in measurements of EWs. 
The sample spectra are characterized by different S/N ratios 
and it is not possible to estimate a typical error in EWs; how- 
ever, errors in the derived abundances due to errors in EWs are 
in a good approximation represented by the standard deviation 
(or rms) around the mean abundance determined from individ- 
ual lines. The rms in principle includes also errors due to un- 
certainties in atomic parameters, but the latter should be mini- 
mized in our analysis, since it is carried out differentially with 
respect to the Sun and our stars have parameters (T e ff, log g) 
close to the solar ones. As already mentioned, rms values for 
[Fe/H] are listed in Tabled while errors in [X/Fe] ratios listed 
in Table @] correspond to the quadratic sum of rms for [Fe/H] 
and rms for [X/H]. 



Errors in [O/H] include in principle uncertainties in the 
measurement of the EW of the [O i] 6300.3 A feature and errors 
due to uncertainties in Ni abundance, which affect the estimate 
of the contribution of the Ni 1 6300.34 blend to the [O i] feature. 
We find however that the latter are much smaller than typical 
errors in [O i] EWs. Our spectra are not characterized by an 
extremely high S/N in the forbidden line region and thus errors 
in EWs are significant. Typical errors are of the order of 0.5- 
1 mA, which reflect into uncertainties in [O/H] between 0.05 
and 0. 1 dex. On the other hand, uncertainties in Ni abundance 
are of the order of 0.05 dex at most and correspond to errors in 
the Ni i EW of the order of 0.15-0.2 mA. 

Internal errors due to uncertainties in stellar parameters 
were estimated by varying each parameter separately, while 
leaving the other two unchanged. We assumed random uncer- 
tainties of ±70 K, +0.15 km/s, and ±0.25 dex in T e ff, £ and 
log g, respectively. We performed different tests and found that 
for all the stars, T e ff variations larger than 70 K would have in- 
troduced significant (i.e, much larger than in the Sun) trends of 
log n(Fe i) vs. EP, while variations in £ larger than 0.15 km/s 
would have resulted in significant trends of log n(Fe) vs. EW. 
This applied not only to Fe, but also to the other elements with 
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several lines covering large dynamical ranges in EP and EWs. 
Finally, differences in log g larger than 0.25 dex would have re- 
sulted in differences between log n(Fe i) and log n(Fe n) larger 
than 0.05 dex. In Table H we list errors in [Fe/H] and [X/Fe] 
ratios due to uncertainties in stellar parameters for the coolest 
(S1239) and one of the warmest (S988) stars in the sample. 

Systematic errors are more difficult to evaluate. Errors in 
the scale of T e ff are likely small, since, as already noted, we did 
not find major trends of inferred Fe abundances as a function 
of EP. By using the T e ff vs. B - V calibration of Alonso et al. 
( 1996J, we would have obtained slightly cooler temperatures, 
with differences ranging between 40 and 60 K, that would have 
given a mean metallicity [Fe/H]=-0.01 (i.e., 0.04 dex below 
our estimate) and almost identical [X/Fe] values. More in gen- 
eral, in order to get an estimate of global systematic errors, 
we analyzed the spectrum of one Hyades member (vB187) ob- 
served with the same UVES set-up as our sample stars (com- 
plete results on the chemical analysis of other Hyades members 
will be reported elsewhere). The results for [Fe/H] and [X/Fe] 
for the Hyades member are listed in Tables^S and|5j our de- 
termination of the metallicity for this star is in perfect agree- 
ment with the classical value for the Hyades (Boesgaard & 
Friel lT990l Paulson et al. l2003l Friel|2006 1; as for the other ele- 
ments, relatively few abundance determinations have been car- 
ried out in the past, in spite of the fact that the Hyades is one of 
the best studied open clusters. In the last line of Tables|4]and[5] 
we list the average [X/Fe] ratios for the Hyades taken from the 
compilation of Friel ( 2006 — Table 1 in that paper). Our [X/Fe] 
ratios for most elements are in good agreement with the values 
from the literature, suggesting that our abundance scale should 
not be affected by major systematic errors. On the other hand, 
for vB 187 we derive an O abundance (and [O/Fe] ratio) consid- 
erably below the average value of Schuler et al. ( 2005 ). We note 
however that the discrepancy is most likely due to differences 
in the EWs of the forbidden line, rather than to the analysis. 
For vB187 we measure an EW of the [O i] 6300.3 A line of 

5 mA to be compared with values of 7-8 mA found by Schuler 
et al. for stars with similar temperatures as vB187. Assuming 
their EWs, we would have inferred an [O/H] similar to their 
values. 

4.3. NLTE effects and 3-D 

As well known, the assumption of LTE may introduce system- 
atic errors and may give origin to spurious abundance trends 
when analyzing stars covering large intervals of effective tem- 
peratures, gravities and metallicities. NLTE effects depend on 
stellar temperature and gravity and should not be a major con- 
cern in the present study, since, as already stressed, we carried 
out a differential analysis with respect to the Sun and most of 
our sample stars are similar to the Sun. In any case, Thevenin 

6 Idiart ( 1999) computed NLTE corrections for Fe i and Fe n 
and showed that NLTE corrections are small for stars with solar 
metallicity or higher. At the metallicity of M 67 the lines used 
for Na determination are marginally affected by NLTE effects 
(e.g., Mashonkina 2000 1: in the temperature and gravity range 
of our sample stars NLTE negative corrections are always be- 



Table 6. Random errors due to uncertainties in stellar parame- 
ters. 



S1239 


:T eff =5477 K, log g=3.6,f = 


1.2 km/s 


A 


AT cff = ±70 


A logg=±0.25 


A£ = ±0.15 




(K) 


dex 


km/s 


[Fe/H]i 


0.06/-0.04 


-0.01/0.02 


-0.05/0.06 


[O/Fe] 


-0.05/0.03 


0.12/-0.14 


0.05/-0.06 


[Na/Fe] 


-0.01/0.00 


-0.03/0.01 


0.03/-0.03 


[Mg/Fe] 


-0.01/0.00 


-0.06/0.05 


0.02/-0.03 


[Al/Fe] 


-0.02/0.00 


0.00/-0.01 


0.04/-0.05 


[Si/Fe] 


-0.05/0.04 


0.03/-0.03 


0.03/-0.04 


[Ca/Fe] 


0.0/-0.02 


-0.04/0.01 


-0.01/0.01 


[Ti/Fe] 


0.02/-0.04 


0.00/-0.01 


0.01/-0.02 


[Cr/Fe] 


0.03/-0.03 


-0.01/0.0 


-0.02/0.02 


[Ni/Fe] 


-0.01/0.01 


0.02/-0.03 


0.01/-0.01 


S988: 


T eft =6151K, log g=4.1,f = 


1.45 km/s 


A 


AT eff = ±70 


A logg=±0.25 


A£ = ±0.15 




(K) 


dex 


km/s 


[Fe/Hi] 


0.05/-0.05 


-0.02/0.02 


-0.03/0.04 


[O/Fe] 


-0.04/0.04 


0.13/-0.14 


0.03/-0.04 


[Na/Fe] 


-0.02/0.01 


-0.02/0.01 


0.01/-0.03 


[Mg/Fe] 


-0.01/0.01 


-0.05/0.05 


0.02/-0.03 


[Al/Fe] 


-0.02/0.02 


0.01/-0.01 


0.02/-0.04 


[Si/Fe] 


-0.02/0.03 


0.02/-0.01 


0.02/-0.03 


[Ca/Fe] 


0.00/0.00 


-0.02/0.02 


-0.01/0.01 


[Ti/Fe] 


0.01/-0.01 


0.01/-0.01 


0.02/-0.02 


[Cr/Fe] 


0.01/-0.02 


0.00/-0.01 


-0.02/0.0 


[Ni/Fe] 


-0.01/0.0 


0.01/-0.02 


0.01/-0.02 



low 0. 1 dex and differential corrections with respect to the Sun 
are of the order of 0.02-0.03 dex at most. NLTE corrections 
are also small for the two Mg lines used in this study (Zhao et 
al. irWSl and the same holds for Al (Baumuller et al.QWSJ. 

Use of time dependent, 3-D hydrodynamical model of the 
solar atmosphere has resulted in the revision of solar abun- 
dances for several elements (Asplund et al. 2005 ). As for NLTE 
effects, use of 1-D models, should not be a major concern for 
our differential analysis. 

4.4. Comparison with other studies 

Our mean [Fe/H] value for M 67 is in very good agree- 
ment with the metallicity determinations of Garcia Lopez et 
al. frgggl [Fe/H] =0.04 + 0.04), Friel & Boesgaard fP^ 
[Fe/H] =0.02+0. 12) and Yongetal. (120051 [Fe/H]= 0.02+0.14), 
while somewhat higher than the values of Tautvaisiene et al. 
(2000] [Fe/H]= -0.03 +0.03) and Shetrone & Sandquist (2000 
[Fe/H]= -0.05), but still consistent with them. The metallicity 
derived by Friel et al. (2002, [Fe/H] =-0.1 5 ± 0.05) from low 
resolution spectroscopy of evolved cluster stars is instead far 
below our estimate. 

To our knowledge, very few studies have been published 
on the abundance pattern of M 67 and significant discrepan- 
cies exist between them. Peterson ( 1992) reports a 0.1 dex a- 
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enhancement based on O, Mg, and Si measurements in two 
cluster giants; on the contrary, Garcia Lopez et al. ( 1988 1, based 
on the analysis of warm unevolved cluster stars, find average 
Ca and Si abundances below solar ([Ca/H]= -0.10 + 0.08 or 
[Ca/Fe]= -0.14, [Si/H] = -0.20 ± 0.06 or [Si/Fe]=-0.24). Both 
the Ca and Si analyses relied on one line only. We have three 
stars in common with the sample of Garcia Lopez et al. {1988 1, 
namely, S988 (F129), S995 (F127), and S2205 (F128); the 
agreement between their and our [Fe/H] values is good, while 
we find higher Ca (significantly higher in the case of F129) and 
much higher Si abundances. 

Tautvaisiene et al. (2000) carried out a detailed chemical 
analysis of 9 clump and red giant branch (RGB) cluster mem- 
bers, finding a rather normal (i.e., close to solar) abundance 
pattern, apart from Na that appeared enhanced. Similar results 
have been obtained more recently by Yong et al. (2005) based 
on three clump cluster members (their stars are in common with 
the sample of Tautvaisiene et al.), although their [X/Fe] ratios 
are slightly higher than those of Tautvaisiene et al. In Table Q 
we compare our average [X/Fe] ratios to those of Tautvaisiene 
et al. for the elements determined in both studies. Specifically, 
we list in Col. 2 the mean abundance ratios from the present 
study and in Cols. 3 and 4 the average ratios from the whole 
sample of Tautvaisiene et al. and those obtained considering 
only the stars observed by them with a resolution R=60,000. 
Focusing on the latter, the table shows that their [X/Fe] ratios 
are in general larger than ours but, considering errors, the abun- 
dance ratios for most elements are consistent with each others. 
Also their [Fe/H] is 0.06 dex below our value, implying that 
the [X/H] ratios are in better agreement. The only exceptions 
are Na, Al, and Si, for which they derive substantially higher 
[X/Fe] (or [X/H]) values than us. 

These discrepancies could be either real, thus indicating in- 
trinsic differences in the chemical composition of unevolved 
and evolved cluster stars and hence chemical processing, or 
due to systematic offsets between Tautvaisiene et al. ( 2000 1 and 
our abundance scale. In order to investigate this point, we re- 
determined Fe, Na, Al, and Si abundances for the three clump 
stars in their sample observed at high resolution, using the same 
code and atomic parameters employed for our sample stars, 
but their stellar parameters. For the analysis, we considered 
only lines in their line list which were also included in ours. 
In Table [8] we compare the [X/Fe] derived by Tautvaisiene et 
al. (|2000) with those obtained with our new analysis of these 
stars: the table shows that, for the three elements and for the 
three stars, we find smaller [X/Fe] ratios than Tautvaisiene et 
al., with differences up to ~ 0.15 dex. Correspondingly, the 
average ratios abundances of Na, Al, and Si to Fe determined 
through our reanalysis are lower than those listed in the last 
column of Tableland are all now consistent with the mean ra- 
tios that we derive for unevolved stars. Sodium remains slightly 
enhanced, but the difference between dwarfs and giants can be 
explained by NLTE effects that are larger for cool giants than 
for warm dwarfs (Mashonkina et al. 2000 1. In other words, our 
analysis suggests that no intrinsic differences are present be- 
tween unevolved and evolved stars in M 67, implying that the 
latter have not undergone significant chemical processing. 



Table 7. Mean abundance ratios from this study and from 
Tautvaisiene et al. (2000). 



Element 


Present 


Tautvaisiene 


Tautvaisiene 


ratio 


study 


all 


R=60,000 


[Fe/H] 


0.03 ±0.03 


-0.03 ±0.03 


-0.02 ±0.03 


[O/Fe] 


0.01 ±0.03 


0.02 ±0.06 


0.05±0.03 


[Na/Fe] 


0.05 ±0.07 


0.19 ±0.06 


0.20 ±0.01 


[Mg/Fe] 


0.00 ±0.02 


0.10 ±0.04 


0.05 ±0.02 


[Al/Fe] 


-0.05 ± 0.04 


0.14 ±0.04 


0.10 ±0.02 


[Si/Fe] 


0.02 ± 0.04 


0.10 ±0.05 


0.08 ±0.02 


[Ca/Fe] 


0.05 ±0.04 


0.04 ±0.05 


0.00 ±0.05 


[Ti/Fe] 


-0.02 ± 0.04 


0.04 ±0.10 


0.03 ±0.02 


[Cr/Fe] 


-0.01 ± 0.04 


0.10 ±0.07 


0.04 ±0.02 


[Ni/Fe] 


-0.02 ± 0.04 


0.04 ±0.04 


0.01 ±0.02 



To our knowledge, this is the first case where an agree- 
ment between abundances in unevolved and evolved stars in 
the same open cluster is found based on a significant sample of 
stars. In particular, the good agreement of Na abundances for 
dwarf and giant stars in the cluster suggests that the slight en- 
hancement with respect to the Sun is probably intrinsic to the 
cluster and not due to the deep mixing in giants as hypothesized 
by Tautvaisiene et al. (|2000) and by Pasquini et al. (2004} for 
IC 4651. 

In summary, the comparison of our abundances for M 67 
with those determined by others shows that systematic offsets 
exist between different studies; most obviously, when investi- 
gating the abundance trends of open clusters as a function of 
age or Galactocentric distances, one should make sure that the 
data are on the same abundance scale. 

5. Discussion and conclusions 

5.1. Star-to-star scatter and Li 

As mentioned in Sect. 1, the investigation of the presence (or 
lack thereof) of a significant dispersion in a and Fe-peak ele- 
ment abundances among cluster stars, that could possibly ex- 
plain the large star-to-star scatter in Li abundances, was the 
main motivation for the present study. 

Figures|2]and|3]together with TablesHI3|5]and|6lclearly in- 
dicate that the standard deviations cr of [Fe/H] and [X/Fe] ratios 
for the whole sample are comparable or even smaller than mea- 
surement uncertainties of individual stars. This on the one hand 
implies that our estimate of internal errors may be somewhat 
conservative; on the other hand, and most important, <x values 
represent upper limits to the dispersion in [Fe/H] and [X/Fe] ra- 
tios, suggesting that we can exclude the presence of a scatter in 
[Fe/H] and [X/Fe] ratios at a level larger than ~ 0.05 dex. In or- 
der to investigate more in detail the relationship between heavy 
element abundances and Li dispersion, sample stars shown in 
Fig.|2]are denoted with different symbols according to whether 
they are Li-rich or Li-poor. Analogously, we show in Fig. 0] 
[X/Fe] vs. [Fe/H] for our sample stars with Li-rich and Li- 
poor objects indicated by different symbols. We define as Li- 
rich/poor stars those lying on the upper/lower envelope of the 
log n(Li) vs. T e tf distribution of M 67 (see Pasquini et al. 119971 
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Table 8. Reanalysis of Tautvaisiene et al. (|2000) data. The subscript "our" indicates the values obtained through our analysis, 
while the subscript "TETI" indicates their original values. 



Ion 


F84 




F141 




F151 




average our 




[X/Fe] TETI 


[X/Fe] our 


[X/Fe] TETI 


[X/Fe] our 


[X/Fe] TETI 


[X/FeU,. 




Nai 


0.21 


0.11 


0.19 


0.12 


0.21 


0.12 


0.12 


All 


0.08 


-0.04 


0.12 


-0.03 


0.09 


-0.01 


-0.03 


Si i 


0.09 


-0.03 


0.06 


-0.03 


0.08 


-0.04 


-0.03 



Jones et al. [T999t . Excluding stars S1034 and S1239 that are 
already on the subgiant branch and have undergone a certain 
amount of post-MS Li dilution (Balachandran 1995 1, three and 
five of the remaining sample stars fall in the Li-poor/rich cat- 
egory (open and filled circles, respectively). The figure shows 
that no systematic difference between the abundance pattern 
of Li-rich and Li-poor stars is present, in the sense that Li- 
poor stars are not systematically more metal-rich/poor and/or 
have larger/lower [X/Fe] ratios than the other stars. In particu- 
lar the two stars S1252 and S1256 (similar temperatures, but a 
factor of almost 10 difference in Li) have very similar [Fe/H] 
and similar (almost identical in some cases) abundances of a- 
elements, including oxygen, as well as of Fe-peak elements. As 
for the other three stars with similar temperatures, but different 
Li abundances (S2205 and S988 both Li-poor and S994 Li- 
rich), they have similar [X/Fe] ratios for most elements — again 
including O — with the exception of Si, Ca, and Cr. [Si/Fe], 
[Ca/Fe] and [Cr/Fe] are somewhat higher in S988 than in the 
other two stars. However, not only the difference is within er- 
rors, but the largest difference in heavy element abundances is 
seen between the two Li-poor stars themselves, with S988 be- 
ing on average more metal rich than S2205. S994 has [X/Fe] 
values intermediate between those of S2205 and S988. 



We conclude that, based on our sample of M 67 members, 
it seems very unlikely that dispersion of a factor about 10 in 
Li abundances measured among otherwise similar cluster stars 
is due to differences larger than ~ 0.05 dex in heavy element 
abundances. We note in particular that no dispersion in O -one 
of the most critical elements for stellar opacities and thus Li 
destruction- is present among our sample stars. 

Our result on the lack of a dispersion in heavy element 
abundances leaves early depletion due to angular momentum 
loss and transport as the most likely explanation for the dis- 
persion (Jones et al. U999t . In this scenario, stars with different 
initial rotational velocities would undergo different amounts of 
Li depletion, resulting in a dispersion in Li at the age of M 67. 
This scenario however encounters two main difficulties: first, 
models including mixing due to angular momentum transport 
predict a correlated Li and Be depletion, which is instead not 
seen in M 67 (Randich et al. 2002 1. Second, since observa- 
tions of young cluster show that they all have similar rotation 
distributions, one would expect that all old clusters should be 
characterized by a spread in Li, which is instead not the case. 
Further studies are most obviously warranted. 



5.2. M 67 in the disk 

In Fig. |5] we compare the average [X/Fe] vs. [Fe/H] pattern of 
M 67 with the sample of field stars thin and thick disks from 
Bensby et al. (EUOH and Bensby et al. d^UOH for O. With the 
caveat that the comparison of M 67 with field stars can be made 
only on qualitative grounds, since abundances are on different 
scales and systematic offsets might be present, the figure in- 
dicates that the average [X/Fe] ratios for M 67 very well fit 
into the general trend of field stars. The c-elements have so- 
lar ratios (or slightly below solar in the case of Ti) and none of 
them is significant over- or under-abundant with respect to field 
stars. [O/Fe] is also solar and very well fits within the distribu- 
tion of field stars. Na, which we find to be somewhat enhanced 
with respect to the Sun, lies on the upper envelope of the dis- 
tribution of field stars, but is still consistent with it. Similarly, 
Al is slightly below the solar ratio and the lower envelope of 
field stars, but, given the uncertainties, consistent with it. In 
other words, our measurements confirm that M 67 has a "nor- 
mal" abundance pattern for the solar neighborhood and it is 
very much similar to the Sun. As a final remark, we wish to 
stress that, given its age and global abundance pattern, M 67 
represents one of the most promising clusters where to look for 
true solar-analogs. 

Acknowledgements. This research has made use of the SIMBAD data 
base, operated at CDS, Strasbourg, France. We thank the anonymous 
referee for the very useful suggestion of including analysis of oxygen 
in the present study. SR and PS acknowledge support by COFIN 2003- 
029437 by the Italian MIUR. 

References 

Allende Prieto, C, Lambert, D.L., and Asplund, M. 2001, ApJ, 556, 
L63 

Alonso, A., Arribas, S., Martinez-Roger, C. 1996, A&AS, 117, 227 
Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 
197 

Asplund, M., Grevesse, N., Sauval, A.J. 2005, in Cosmic Abundances 
as Records of Stellar Evolution and Nucleosynthesis, Thomas G. 
Barnes III and Frank N. Bash (eds), ASP Conf. Ser. 336, p. 25 

Balachandran, S. 1995, ApJ, 446, 203 

Baumuller, D., Butler, K, and Geheren, T. 1998, A&A, 338, 637 
Bensby, T, Feltzing, S., and Lundstrom, I. 2003, A&A, 410, 527 
Bensby, T, Feltzing, S., and Lundstrom, I. 2004, A&A, 415, 155 
Biemont, E., Baudoux, M., Kurucz, R.L., Ansbacher, W., and 

Pinnington, E.H. 1991, A&A, 249, 539 
Boesgaard, A.M., and Friel, E.D 1990, ApJ, 351, 467 
Chen, Y.Q., Nissen, P.E., Zhao, H.W., Zhang, H.W., and Benoni, T. 

2000, A&AS, 141,491 



S. Randich et aL: Chemical abundances in M 67 



9 



Clementini, G., Carretta, E., Merighi, R. et al. 1995, AJ, 110, 2319 
Dekker, H., D' Odorico, S., Kaufer, A., Delabre, B., Kotzlowski, H. 

2002, SPIE 4008, 534 
Friel, E.D. 2006, in Chemical Abundances and Mixing in Stars in the 

Milky Way and its Satellites, L. Pasquini and S. Randich (eds), 

Springer- Verlag, in press 
Friel, E.D. & Boesgaard, A.M. 1992, ApJ, 387, 170 
Friel, E.D., Janes, K.A., Tavarez, M., Katsanis, R., Lotz, J., et al. 2002, 

AJ, 124, 2693 
Fulbright, J.P. 2000, ApJ, 120, 1841 

Garcia Lopez, R.J., Rebolo, R., and Beckmann, J.E. 1988, PASP, 100, 
1489 

Grevesse, N., & Sauval, A.J. 1998, SSRv, 85, 161 

Jeffries, R. 2006, in Chemical Abundances and Mixing in Stars in the 

Milky Way and its Satellites, L. Pasquini and S. Randich (eds), 

Springer- Verlag, in press 
Johansson, S., Litzen, U., Lundberg, H., Zhang, Z. 2003, ApJ, 584, 

L107 

Johnson, H.L. 1966, ARA&A, 4, 193 

Jones, B.F, Fisher, D., and Soderblom, D.R. 1999, AJ, 117, 330 
Kupka, F, Piskunov, N.E., Ryabchikova, T.A, Stempels, H.C., and 

Weiss, W.W. 1999, A&AS 138, 119 
Kurucz, R.L. 1993, ATLAS9 Stellar Atmosphere Programs (Kurucz 

CD-ROM No. 13) 
Lambert, D.L., and Warner, B. 1968, MNRAS, 138, 181 
Lambert, DJ. 1978, MNRAS, 182, 249 

Mashonkina, L.I., Shimanskii, V.V, and Sakhibullin, N.A. 2000, 

Astronomy Reports, 44, 790 
Montgomery, K. A., Marschall, L. A., & Janes, K. A. 1993, AJ 106, 

18 

Nissen, RE. 1981, A&A, 97, 145 

Pace, G, Pasquini, L., and Ortolani, S. 2004, A&A, 401, 997 
Pallavicini, R., Cerruti-Sola, M., Duncan, D.K. 1987, A&A, 174, 116 
Pallavicini, R., Spano, P., Prisinzano, L., Randich, S., and Sestito, P. 
2006, in Chemical Abundances and Mixing in Stars in the Milky 
Way and its Satellites, L. Pasquini and S. Randich (eds), Springer- 
Verlag, in press 

Pasquini, L., Liu, Q., and Pallavicini, R. 1994, A&A, 287, 191 
Pasquini, L., Randich, S., and Pallavicini, R. 1997, A&A, 325, 535 
Pasquini, L., Randich, S., Zoccali, M. et al. 2004, A&A, 424, 951 
Paulson, D.B., Sneden, C, and Cochran, W.D. 2003, AJ, 125, 318 
Peterson, R.C. 1992, AAS, 181, 2307 
Piau, L., and Turck-Chieze, S. 2002, ApJ, 566, 419 
Piau, L., Randich, S., and Palla, F. 2003, A&A, 408, 1037 
Piskunov, N.E., Kupka, F, Ryabchikova, T.A, Weiss, W.W., and 

Jeffery, C.S. 1995, A&AS 112, 525 
Randich, S. 2006, in Chemical Abundances and Mixing in Stars in the 

Milky Way and its Satellites, L. Pasquini and S. Randich (eds), 

Springer- Verlag, in press 
Randich, S., Pasquini, L., Pallavicini, R. 2000, A&A, 356, L25 
Randich, S., Primas, F, Pasquini, L., and Pallavicini, R. 2002, A&A, 

387, 222 

Randich, S., Sestito, P., Pallavicini, R. 2003, A&A, 399, 133 
Ryabchikova, T.A., Piskunov, N.E., Stempels, H.C., Kupka, F, and 

Weiss, W.W. 1999, Physica Scripta T83, 162-173 (VALD-2) 
Sanders, W.L. 1977, A&A Suppl., 27, 89 
Sandquist, E.L. 2004, MNRAS, 347, 101 

Schuler, S.C., Hatzes, A.P., King, J.R., Kurster, M., and The, L. 2005, 

ApJ, in press l astro-ph/05 10721 1 
Shetrone, M.D., & Sandquist, E.L. 2000, AJ, 120, 1913 
Sestito, P., Randich, S., and Pallavicini, R. 2004, A&A, 426, 809 
Sestito, P., degli Innocenti, S., Prada-Moroni, P., and Randich, S. 2006, 

A&A, submitted 



Sneden, C. 1973, ApJ, 184, 839 

Spite, F, Spite, M., Peterson, R. C, Chaffee, F. H., Jr. 1987, A&A, 
171, L8 

Soderblom, D.R., Jones, B.F, Balachandran, S., et al. 1993, AJ, 106, 
1059 

Swenson, F.J., and Faulkner, J. 1992, ApJ, 395, 654 

Tautvaisiene, G., Edvardsson, B., Tuominen, I., and Ilyin, I. 2000, 

A&A 360, 499 
Thevenin, F, and Idiart, TP. 1999, ApJ, 521, 753 
Unsold., A. 1955, Physik der Sternatmospharen. Springer- Verlag, 

Berlin 

Yong, D., Carney, B.W, Teixeira de Almeida, M.L. 2005, AJ, 130, 
597 

Zhao, G., Butler, K, Gehren, T. 1998, A&A, 333, 219 



10 



S. Randich et al.: Chemical abundances in M 67 




B-V 



Fig. 1. Color-magnitude diagram of M 67. The sample stars are indicated as filled circles. CCD photometry was taken from 
Montgomery et al. dl993> . Only stars brighter than V=16 are plotted in the figure. 
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Fig. 2. [Fe/H] as a function of effective temperature for the sample stars. Filled and open circles represent main sequence Li-rich 
and Li-poor stars, while open squares denote the two slightly evolved stars. Error bars shown in the figure correspond to rms 
values quoted in Tabled Internal errors due to uncertainties in stellar parameters are not shown here. The solid and dotted lines 
represent the average [Fe/H] value +1 standard deviation. 
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Fig. 3. [X/Fe] as a function of effective temperature for the sample stars. Error bars correspond to the quadratic sum of the rms 
of [Fe/H] (or errors due to uncertainties in measured EWs for O, whose abundance has been determined from one line only) and 
rms of [X/H]. As in Fig.|3] in each panel the solid and dotted lines represent the mean value +1 standard deviation. 
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Fig. 4. [X/Fe] vs. [Fe/H] for Li-poor (open circles) and Li-rich (filled circles) stars. Error bars in [Fe/H] and [X/Fe] are the same 
as in Figs.|3]and[3] Horizontal lines are the same as in Fig. [5] The two evolved stars are not plotted in the figure (see text). 
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Fig. 5. [X/Fe] vs. [Fe/H] for field stars from Bensby et al. ( 20031 120041 — open circles) and the average values from this study 
(filled symbols). 
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Table 2. Line list, adopted log gf values, and solar EWs. (1): VALD; (2): Lambert & Walner fT9fj8l ; (3): Chen et al. J2000l : (4): Bensby et al. 
J20031 : (5): Clementini et al. Jl995l : (6): Fulbright <2000t : (7): Biemont et al (19911 ; (*): inverse solar analysis. 



Element 


A ( A) 


log gf 


ref. 


EW (mA) 


Oi 


6300.633 


-9.717 


(1) 


5.5 


Nai 


5682.633 


-0.700 


(1) 


107.0 


Nai 


5688.220 


-0.400 


(2) 


130.0 


Nai 


6154.226 


-1.570 


(2) 


35.0 


Nai 


6160.747 


-1.270 


(2) 


52.9 


Mgi 


5528.405 


-0.620 


(1) 


233.0 


Mgi 


5711.090 


-1.724 


(3) 


105.0 


All 


5557.070 


-2.160 


(*) 


11.5 


All 


6696.023 


-1.540 


(*) 


39.3 


All 


6698.673 


-1.890 


(*) 


20.1 


Si i 


5701.104 


-2.050 


(1) 


37.2 


Si i 


5948.541 


-1.230 


(1) 


81.5 


Si i 


6091.919 


-1.400 


(1) 


30.6 


Si i 


6125.021 


-1.550 


(*) 


30.9 


Si i 


6142.483 


-1.480 


(6) 


33.2 


Si i 


6145.016 


-1.440 


(5) 


37.5 


Si i 


6414.980 


-1.100 


(1) 


46.0 


Si i 


6518.733 


-1.500 


(1) 


20.5 


Si i 


6555.463 


-1.000 


(1) 


43.6 


Cai 


5512.980 


-0.480 


(*) 


87.0 


Cai 


5581.965 


-0.671 


(3) 


94.3 


Cai 


5601.277 


-0.523 


(3) 


105.5 


Cai 


5867.562 


-1.610 


(*) 


24.4 


Cai 


6102.723 


-0.862 


(1) 


126.5 


Cai 


6122.217 


-0.386 


(1) 


171.2 


Cai 


6161.297 


-1.293 


(1) 


60.5 


Cai 


6166.439 


-1.156 


(1) 


68.6 


Cai 


6169.042 


-0.804 


(1) 


87.0 


Cai 


6169.563 


-0.527 


(1) 


107.6 


Cai 


6455.598 


-1.400 


(*) 


56.6 


Cai 


6499.650 


-0.818 


(3) 


87.0 


Tii 


4805.415 


0.150 


(1) 


32.6 


Tii 


4820.411 


-0.441 


(1) 


44.9 


Tii 


4885.079 


0.358 


(1) 


64.0 


Tii 


4913.614 


0.160 


(1) 


53.7 


Tii 


5016.161 


-0.574 


(1) 


67.5 


Tii 


5219.702 


-2.292 


(1) 


28.1 


Tii 


5866.451 


-0.840 


(1) 


48.1 


Tii 


5953.160 


-0.329 


(1) 


33.6 


Tii 


5965.828 


-0.409 


(1) 


29.9 


Tii 


6258.102 


-0.431 


(3) 


50.9 


Tii 


6261.098 


-0.479 


(1) 


49.8 


Cri 


4936.335 


-0.340 


(1) 


42.0 


Cri 


5247.566 


-1.640 


(1) 


81.2 


Cri 


5296.691 


-1.400 


(1) 


92.2 


Cri 


5300.744 


-2.120 


(1) 


59.8 


Cri 


5329.142 


-0.064 


(1) 


66.0 


Cri 


5348.312 


-1.290 


(1) 


96.9 


Fei 


4835.868 


-1.500 


(1) 


46.6 


Fei 


4875.878 


-2.020 


(1) 


55.6 


Fei 


4907.732 


-1.840 


(1) 


58.7 


Fei 


4999.113 


-1.740 


(1) 


31.2 


Fei 


5036.922 


-3.068 


(1) 


22.2 



Table 2. continued. 
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Element 


o 

A (A) 


log gf 


ref. 


EW (mA) 


Fe i 


5044.211 


-2.038 


(1) 


75.1 


Fe i 


5067.150 


-0.970 


(1) 


70.0 


Fe i 


5141.739 


-2.190 


(1) 


88.8 


Fe i 


5162.273 


0.020 


(1) 


138.1 


Fe i 


5217.389 


-1.070 


(1) 


110.9 


Fei 


5228.377 


-1.290 


(1) 


52.1 


Fei 


5285.129 


-1.640 


(1) 


26.7 


Fei 


5293.959 


-1.870 


(1) 


29.0 


Fe i 


5373.709 


-0.860 


(1) 


63.0 


Fe i 


5386.334 


-1.770 


(1) 


32.0 


Fei 


5397.618 


-2.480 


(1) 


23.1 


Fe i 


5472.709 


-1.495 


(1) 


42.8 


Fei 


5522.447 


-1.550 


(1) 


41.4 


Fei 


5539.280 


-2.660 


(1) 


16.2 


Fe i 


5543.150 


-1.570 


(1) 


63.3 


Fei 


5543.936 


-1.140 


(1) 


60.2 


Fei 


5546.992 


-1.910 


(1) 


23.9 


Fei 


5584.765 


-2.320 


(1) 


34.8 


Fei 


5636.696 


-2.610 


(1) 


19.1 


Fe i 


5638.262 


-0.870 


(1) 


73.8 


Fe i 


5662.516 


-0.573 


(1) 


93.4 


Fei 


5691.497 


-1.520 


(1) 


39.6 


Fei 


5701.545 


-2.216 


(1) 


84.8 


Fei 


5862.353 


-0.058 


(1) 


86.4 


Fe i 


5916.247 


-2.994 


(1) 


54.2 


Fei 


5930.180 


-0.230 


(1) 


88.8 


Fei 


5934.655 


-1.170 


(1) 


72.6 


Fei 


5956.694 


-4.605 


(1) 


54.1 


Fei 


5976.775 


-1.310 


(1) 


66.7 


Fei 


5984.814 


-0.343 


(1) 


82.7 


Fe i 


5987.066 


-0.556 


(1) 


66.8 


Fe i 


6003.012 


-1.120 


(1) 


77.8 


Fei 


6024.058 


-0.120 


(1) 


108.7 


Fei 


6056.005 


-0.460 


(1) 


72.6 


Fei 


6078.491 


-0.424 


(1) 


76.0 


Fe i 


6136.995 


-2.950 


(1) 


68.6 


Fe i 


6157.728 


-1.260 


(1) 


61.8 


Fei 


6187.990 


-1.720 


(1) 


46.2 


Fei 


6200.313 


-2.437 


(1) 


73.5 


Fei 


6315.811 


-1.710 


(1) 


40.3 


Fei 


6322.685 


-2.426 


(1) 


75.0 


Fe i 


6336.824 


-0.856 


(1) 


105.5 


Fei 


6344.149 


-2.923 


(1) 


61.0 


Fei 


6469.193 


-0.770 


(1) 


54.8 


Fei 


6495.742 


-0.940 


(1) 


42.2 


Fe i 


6498.939 


-4.699 




46.6 


Fei 


6574.228 


-5.023 




29.4 


Fei 


6609.110 


-2.692 




66.2 


Fei 


6703.567 


-3.160 




38.1 


Fei 


6733.151 


-1.580 




23.1 


Fei 


6750.153 


-2.621 




74.2 


Fei 


6806.845 


-3.210 




33.7 


Fe ii 


5264.812 


-3.120 




47.1 


Fe ii 


5325.553 


-3.222 


(7) 


41.1 


Fe ii 


5414.073 


-3.750 


(7) 


23.7 



Table 2. continued. 
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Element 


o 

A ( A) 


log gf 


ref. 


EW (mA) 


Fe ii 


5425.257 


-3.372 


(7) 


37.8 


Fe ii 


5991.376 


-3.557 


(7) 


31.4 


Fe n 


6084.111 


-3.808 


(7) 


21.5 


Fe n 


6149.258 


-2.724 


(7) 


34.8 


Fe ii 


6247.557 


-2.329 


(7) 


53.6 


Fe ii 


6432.680 


-3.708 


(7) 


39.4 


Fe ii 


6456.383 


-2.075 


(7) 


64.4 


Fe ii 


6516.080 


-3.450 


(7) 


49.3 


Ni i 


4806.984 


-0.640 


(1) 


59.2 


Ni i 


4852.547 


-1.070 


(1) 


42.7 


Ni i 


4904.407 


-0.170 


(1) 


87.9 


Ni i 


4913.968 


-0.630 


(1) 


58.7 


Ni i 


4946.029 


-1.290 


(1) 


23.5 


Nil 


5003.734 


-3.130 


(4) 


32.8 


Nil 


5010.934 


-0.870 


(1) 


52.2 


Nil 


5032.723 


-1.270 


(1) 


20.6 


Nil 


5082.339 


-0.590 


(*) 


62.8 


Nil 


5155.125 


-0.650 


(1) 


47.9 


Nil 


5435.855 


-2.590 


(1) 


45.8 


Ni i 


5462.485 


-0.930 


(1) 


39.9 


Ni i 


5589.357 


-1.140 


(1) 


27.1 


Nil 


5593.733 


-0.840 


(1) 


41.6 


Ni i 


C £. Hi C Oil 

5625.312 


-0.700 




38.6 


Nil 


5641.880 


-1.070 




24.1 


Nil 


5682.198 


-0.499 




50.1 


Nil 


6111.066 


-0.870 




33.4 


Nil 


6175.360 


-0.559 




48.4 


Nil 


6186.709 


-0.960 




29.0 


Nil 


6191.171 


-2.353 




74.9 


Nil 


6223.981 


-0.910 




26.2 


Nil 


6378.247 


-0.830 




33.3 


Nil 


6586.308 


-2.810 




41.0 



